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^ . Abstract. We explore the fate of the universe given the possibility that the density 

^ il associated with 'dark energy' may decay slowly with time. Decaying dark energy 

' is modeled by a homogeneous scalar field which couples minimally to gravity and 

, whose potential has at least one local quadratic maximum. Dark energy decays as 

^SJ ' the scalar field rolls down its potential, consequently the current acceleration epoch is 

a transient. We examine two models of decaying dark energy. In the first, the dark 

I energy potential is modeled by an analytical form which is generic close to the potential 

' maximum. The second potential is the cosine, which can become negative as the field 

. evolves, ensuring that a spatially flat universe collapses in the future. We examine the 

I feasibility of both models using observations of high redshift type la supernovae. A 

■ maximum likelihood analysis is used to find allowed regions in the {m, 0o} plane (m is 

fT^ , the tachyon mass modulus and (f>Q the initial scalar field value; m ^ Hq and 00 ^ Mp 

' by order of magnitude). For the first model, the time for the potential to drop to half 

. its maximum value is larger than ~ 8 Gyrs. In the case of the cosine potential, the 

Oh| time left until the universe collapses is always greater than ~ 18 Gyrs (both estimates 

Q ■ are presented for l^om = 0.3, m/Ho ^ 1, Hq ~ 70 km/sec/Mpc, and at the 95.4% 

^ , confidence level). 



, 1. Introduction 



Observations of the luminosity distance to distant type la supernovae Pfl Ej supported 
by the discovery of CMB angular temperature fluctuations on degree scales jS] and 
measurements of the power spectrum of galaxy clustering |3] convincingly show that our 
Universe is (approximately) spatially fiat, with ~ 30% of its critical energy density in 
non-relativistic matter (cold dark matter (CDM) and baryons) and the remaining ~ 70% 
in a smooth component having a large negative pressure ('dark energy'). Dark energy is 
clearly the most abundant form of matter in the Universe (in terms of the effective energy 
density), yet both its nature and its cosmological origin remain enigmatic at present. It 
is clear that the observational properties of dark energy (DE) are tantalizingly close to 
those of a cosmological constant 7', yet dark energy need not be / exactly (see reviews 

0ini). 

An intriguing possibility, carrying far reaching consequences both for the current 
Universe and its ultimate fate, is that the dark energy may be decaying. The possibility 
that dark energy could be unstable is in fact suggested by the remarkable qualitative 
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analogy between the presence of dark energy today and the properties of a different type 
of 'dark energy' - the inflaton field - postulated in the inflationary scenario of the early 
Universe. This analogy works in two ways. On one hand, the fact that a form of matter 
having a large negative pressure dominates the universe today makes it not unnatural 
that a similar form of matter (having w < 0) could have dominated the universe in the 
distant past. On the other hand, since dark energy in the early Universe (the inflaton) 
was unstable and decayed aeons ago, one might be tempted to ask whether the nature 
of dark energy observed today will be any different. In this paper we address in detail 
the possibility that dark energy, like the inflaton which preceded it, may be decaying. 

Decaying dark energy (DDE) leads to numerous interesting possibilities including 
the fact that the current epoch of cosmic acceleration could be a transient which ends 
after the dark energy density has dropped to sufficiently small values. Such a universe 
will clearly be very different from the standard ACDM cosmology. For instance, a DDE 
universe may not possess horizons which are characteristic of ACDM as well as tracker- 
driven quintessence models (see |Zj for a brief review). (Interesting implications for the 
future of our universe also arise in certain braneworld models, in which the current 
accelerating regime is a transient between two matter dominated epochs 1^1 llOj.) 

Within the context of quintessence models another interesting (though more 
speculative) possibility is provided by potentials which are not constrained to be positive 
but which can become negative for certain values of ^T]. These potentials are either 
bounded from below, in which case V{(f)) has one or more minima at which V{(j)i) < 0. 
(A good example is furnished by the cosine potential which will be examined in detail 
later in this paper.) Alternatively the potential becomes unbounded, in which case 
y(0) ^ —oo, for some values of (p. Under the influence of a negative potential, DDE 
steadily decreases and the universe expands at an increasingly slower rate. Finally a 
stage is reached when the negative DDE density exactly cancels the positive density of 
dark matter (plus the curvature term, if present). At this point of time the expansion of 
the universe stops {H = 0) and thereafter the universe begins to collapse {H < 0). This 
paper examines two general classes of potentials describing decaying dark energy in the 
light of recent high redshift supernova observations. In the first class V{(f)) = Vq — |m^0^, 
and we evolve the field equations into the future to determine the time taken for the 
potential to drop to half its maximum value. The second class of models involves the 
cosine potential. In both cases we use the agreement between our model and high 
redshift supernova data in the past, to extrapolate and predict the behaviour of our 
model universe in the future. As a result, we find the range of model parameters allowed 
by the supernovae data and determine: (i) the minimum time for the universe to stop 
accelerating, and, (ii) the minimum time for the universe to collapse in DDE models 
with negative potentials. 
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2. Decaying Dark Energy 

The decaying dark energy (DDE) model which we consider is analogous to successful 
inflationary models in that it undergoes quasi-homogeneous decay We model DDE by 
a scalar field (j) ('quintessence') which couples minimally to gravity and has one or more 
maxima in its potential. Near any one such local maximum (which for convenience we 
assume is at = 0) the field can be generically modeled by the potential 



V{4>) = Vq — , where = ' 



(1) 

(t>=0 



The modulus of the tachyonic rest- mass m will be assumed ~ Hq where Hq is the Hubble 

constant (the slow-roll condition, m <^ Hq, would result in the absence of any noticeable 
observational effects for ultra-light scalars, making such fields virtually indistinguishable 
from ACDM, at least within the accuracy of present data). It is known from observations 
that physical properties of dark energy are close to those of a cosmological constant at 
present. This means that the variation in must be on time scales of order Hq^ and not 
much faster. (In the opposite case the equation of state w would vary much too rapidly.) 
The equation of motion then implies that m ~ ifo which suggests that all characteristic 
time-scales in our problem will be of order Hq^ or greater. Therefore, in this model the 
initial value 0o oi the scalar field (p and its change - during the entire process of dark 
energy decay - are both expected to be of the order of G~^. (This follows quite simply 
from the Einstein equation m^^^ ~ p^^.^^ ~ Hq/G.) However, since cosmology is now a 
precision science, we can hope to obtain much better quantitative results, beyond these 
simple qualitative estimates, using current observational data. This is precisely what is 
done in this paper. 

We should at this stage mention that, in principle, decay mechanisms other than 
those considered in this paper are possible. As an example one should mention bubble- 
like decay via quantum tunneling in the case when V" > 0, or strongly inhomogeneous 
classical decay in the case of a tachyonic mass large compared to Hq, however these 
possibilities must be confronted with observational data other than the high-z supernova 
data which we wish to consider, and will therefore not be discussed any further by us 
in this paper. It may also be appropriate to mention that our model of decaying dark 
energy does not introduce any direct non-gravitational coupling of dark energy to non- 
relativistic matter. Significant coupling of dark energy to dark matter may lead to large 
energy transfer from dark energy to dark matter in the process of DDE decay. However 
this model of DDE, though potentially interesting, lies somewhat beyond the scope of 
the present paper and we do not consider it any further here. In our case, what we call 
the DE decay is actually the transition of DE from the potential-dominated regime to 
a regime in which the potential and kinetic energies of the dark energy field (j) become 
comparable. | 

I Wc also do not consider fields which couple non-minimally to gravity. Such a coupling, if small, would 
result in an insignificant change in the eflFective mass parameter m for time intervals of the order of 
Hq^ which we consider in this paper. Large couplings would lead to a value of the effective Brans-Dicke 
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The field equations governing the behaviour of the scalar field as it rolls down its 
potential in a spatially flat universe are 

4> +3H,p + V\(f)) = , 

H'=^-^ {Pr^ + p,) , (2) 

where V{(j)) is given by (Q) and 

P0 = V^(0) + ■ (3) 

At early times {Hot -C 1), the universe was matter-dominated. Assuming p^ ^ p</, in 
(0) we obtain the following exact solution for (f){t) during the matter dominated regime: 

, , , , sinh mt , , , 

t = 00 — . 4 

mt 

Eq (0)) determines the initial conditions for and 0. At very early times {mt <^ 1) the 
very large damping experienced by the scalar field ensures that it remains close to its 
initial value at 0o. One should note that this initial value must be sufficiently small, 



00 <, yS/SvrG, in order that DE have sufficiently large negative pressure during the 
present epoch. This is clearly demonstrated in the next section when we compare this 
model of dark energy with supernova observations. 

Since the point = 0, = is a saddle point for homogeneous solutions of the 
scalar field equation with the potential (0), there exists a set of non-zero measure of 
generic solutions which remain sufficiently close to our solution up until the present 
epoch. (It is sufficient for this purpose that a large initial kinetic energy - if at all it 
existed - is redshifted by ^ ~ 3 so that the field settles on the trajectory described by 
(0)) by that redshift. Since the kinetic energy of the scalar field decays as oc during 
the regime when 0^ ^ m^0^ this requirement is easily satisfied in practice.) Of course, 
one may with good reason ask as to whether our set of initial conditions on {0, 0} is not 
too small. The answer to this question depends entirely on the (unknown) behaviour of 
present dark energy at large temperatures and curvatures where the form of the potential 
need not be valid. For this reason, we will not discuss it any further here. However 
it should be noted that, generally speaking, the probability of such initial conditions 
will be strongly enhanced if there are many maxima in the DE potential. This occurs, 
for instance, in the DE model with V(0) = Vq cos(0//) +Vi [IS CHI CSl E! which will 
be discussed by us later in this paper (in the inflationary context, such a model is called 
'natural inflation'). 

In the next section, we numerically integrate the system of equations Q with the 
potential (HJ, and determine the scale factor a{t) and the Hubble parameter H{z) = a/a 
as functions of the initial field displacement 0o for specific values of m and Qom- Our 

parameter to which violates the lower bounds on this parameter set by Solar system tests of Einstein 
gravity. 
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model is confronted against high redshift type la supernova observations through its 
luminosity distance 

which is used to place constraints on the free parameters of the model, namely m and 
00 {Vq is not a free parameter since it is uniquely determined by m, (pQ, Qq^ and the 
current value of the Hubble parameter - Hq). 

Having determined the permitted range of parameters values, we shall proceed to 
determine the future evolution of our model universe. As emphasized in ^3], reliable 
future predictions are only possible for finite intervals of time. Previous predictions 
(of the kind 'our Universe will keep expanding and not encounter a singularity') had a 
depth of approximately 20 Gyrs jHj (see [El UHl CHI for similar estimates and [Oj for the 
braneworld context). It is entirely reasonable to use the form ^ up to the point when 
m^0^ = Vq. This occurs when the potential has declined to half its maximum value. In 
the next section we shall determine the minimal period of time for the potential to reach 
its half-way mark (this time could be referred to as the 'half-decay time of dark energy'). 
Note that this estimate is rather robust since it depends upon the fairly general form of 
the potential near its maximum value given by and not on any other details of dark 
energy. 

In Sec. 3, we also consider a potential which can become negative. The form of 
this potential is assumed to be fHi IT^ ITEj : 



\/ = \/ocos^, f = ^ (6) 

where the value of / is chosen in such a way that the potential (jHI) coincides with ^ 
at small 0. Of course, the hypothesis that V may acquire negative values is speculative 
since it does not follow from any of the current observational data. However, such 
potentials often arise in supergravity and M-theory models (see, e.g., jT3 HZl d]). For 
this potential, we calculate: (i) the minimal extent of the current acceleration epoch, and 
(ii) the time elapsed before the universe collapses. (We should draw the readers attention 
to the fact that the collapse of the universe is a generic property of flat cosmological 
models with negative potentials [H].) We will not however consider the subsequent 
period of contraction in any detail since its exact duration strongly depends both upon 
the behaviour of V in the region ^ < and upon the properties of dark matter at 
high energies. We end with section 4 which contains a summary of our results and a 
discussion. 



3. Methodology and Results 

We constrain the parameter space of our cosmology by requiring that our dark energy 
model provides a good fit to type la supernova data. For this purpose we use the 54 
SNe la from the primary 'fit C of the Supernova Cosmology Project, which includes 16 
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low redshift Calan-Tololo SNe PQ . Fit C is a subsample of a total of 60 SNe of which six 
are excluded as outhers: two low redshift SNe due to suspected reddening and four high 
redshift SNe which are excluded due to atypical light curves. The measured quantity 
in this data, the bolometric magnitude m^, is related to the luminosity distance and 
therefore the cosmological parameters by the following equation 



where Dl = H^di is the Hubble-parameter-free luminosity distance and A4 = Mb + 
25 — 5 log^o-^o is the Hubble-parameter-free absolute magnitude. We shall assume 
that the SNe measurements come with uncorrelated Gaussian errors in which case the 
likelihood function is given by the chi-squared distribution with N—n degrees of freedom: 
C oc exp (— (In our case, = 54 and n = 3). 
The x^-statistic is defined as 



where is the effective B-band magnitude of the i-th supernova obtained after 
correcting the observed magnitude at redshift z for the supernova width-luminosity 
relation, am, is the error in magnitude at redshift z, and m{zi) is the apparent magnitude 
of the i-th supernova in our cosmological model. 

For a given value of f2om, the parameters to be estimated in our cosmology are Ai, 
m and 0o- For practical purposes, the quantity is a statistical nuisance parameter, 
and we marginalize over it (assuming a uniform prior) when we determine m, 0o- We 
perform the minimization while evolving the scalar field according to the equations of 
motion For this purpose, we shall find it convenient to measure the parameter 
m in units of the quantity wSttCVo/S and the parameter 0o in units of the reduced 



Planck mass Mp = yS/STrG. (In these units we constrain the value of m using the 
prior m > 0.1. This condition is set to minimize numerical uncertainties and to avoid 
models which are virtually indistinguishable from ACDM). 

In Figure ^ we show the confidence levels in the m — 0o plane for the first potential 
described by Eq (^. We consider two cosmological models in which the current value 
of f2om is 0.3 and 0.4 respectively. (The results for the two cases differ marginally.) We 
see that an enormous region of parameter space which corresponds to large values of 
m & 00 is disallowed in both models. The reason for this is simple, for large values 
m ^ 1, 00 ^ 1; the scalar field potential decreases much too rapidly and the matter 
density never gets to reach its present value of fiom = 0.3 (left panel) or Qq^ = 0.4 
(left panel). Also, from the confidence levels we see that at higher values of m {(po), the 
model provides a good fit to SNe data only if the corresponding value of 0o (m) is very 
small. (It should be noted that for asymptotically small values of either m or 0o the 
DDE model becomes virtually indistinguishable from ZCDM.) From the value of the 
at best-fit we see that higher values of f2om seem to be slightly favoured by models of 
DDE. For Qom = 0.3, xL = 1-053 for the best-fit model with m = 0.73, 0o = 0.18, while 



mB = M + 5 log^QDiiz; Qq^, m, 0o) , 



(7) 




(8) 





Figure 1. Confidence levels at 68.3% (light grey inner contour) 95.4% (medium 
grey contour) and 99.73% (dark grey outer contour) are shown in the m-(j)Q 
plane for the potential F(</>) =Vq — ^m^cjp'. In the left panel, the present value 
of the matter density is Oom = 0.3, and in the right panel it is = 0.4. The 
mass m is measured in units of \/87rGVo/3, and the initial field value, (f)^, is in 
units of the reduced Planck mass Mp = ^/SjSnG. The best-fit point in each 
plot is marked by a star. The psr degree of freedom at the best-fit is 1.053 
for r^om = 0.3 and 1.049 for J7om = 0.4. For both figures, in the region to the 
right of the thick solid line, parameter values are such that the matter density 
never reaches its present value, hence this region is disallowed by observations. 



/CDM (for identical fiom) has xL = 1-054. For fiom = 0.4, xL = 1-049 for the best-fit 
model m = 0.81, 0o = 0.19, while Xdof(^CDM) = 1.058 (for identical fiom)- For lower 
values of f2om, ^-g- ^om = 0.2, /CDM is marginally favoured over the best-fit DDE. 

In Figure El we examine this potential in greater detail by restricting ourselves to 
a smaller (and more probable) region in parameter space. In this region we show the 
lines of constant AT1/2, where AT1/2 is the time taken, as measured from the present 
epoch, for the potential to drop to half its maximum value: V[(t>itQ + AT1/2)] = Vq/2. 
We see that, at the 95.4% confidence level, the minimum time at which this happens is 
ATi/a ~ 8 Gyrs (if Hq ~ 70 km/sec/Mpc). 

In Figure El we show the confidence levels in the m — <^o plane for the second 
potential, described by Eq The results in this case are qualitatively similar to those 
for the first potential. Here too, a large region of parameter space corresponding to 
higher values of m and (po is disallowed because for these parameter values the matter 
density never reaches its present day value. Also, at higher values of m (^o); the model 
provides a good fit to SNe data only if the complementary parameter 0o {^) is very 
small. Again, for asymptotically small values of either m or 0o this decaying dark 




Figure 2. A magnified part of tlie Fig. ^ witli (dashed) lines of constant 
added. Ari/2 is the time, measured from the present epoch, to when 
the DDE potential has dropped to half its maximum value: V{(p) = Vo/2. 
The values of ATi^2 the dashed curves (from top to bottom) are listed 
in Table [T] (from left to right). For both J7om = 0.3 and Oom = 0.4, the 
minimum time elapsed before the potential drops to half its maximum value is 
ATi/2 - ~ 8 Gyrs(i/o = 70 km/s/Mpc) at the 95.4% confidence level. 

In the region to the right of the thick solid curve, parameter values are such 
that the matter density never reaches its present value. This region is therefore 
disallowed by observations. 



Table 1. Time taken for the potential to drop to half its maximum value. 





ATi/2 (Gyr) 


0.3 


8.3 


42.0 139.7 


712.8 


0.4 


8.3 


34.9 153.6 


712.8 



energy model becomes virtually indistinguishable from ZCDM. For fiom = 0.3, we have 
Xdof = 1-050 at the best-fit point of m = 0.74, 0o = 0.23, and for fiom = 0.4, xL = 1-047 
for the best-fit model having m = 0.79, 0o = 0.24. As in the previous analysis, ZCDM is 
marginally preferred over the best-fit DDE for Qq^ < 0.3. 

In Figure m we examine the results for this potential more closely by focusing on a 
smaller region in parameter space. In this region we show the lines of constant ATend, 
which is the time left until the universe ceases to accelerate. (In terms of the deceleration 
parameter, q{to + ATend) = 0.) We see that, at the 95.4% confidence level, the universe 
will continue accelerating for at least ~ 10 Gyrs. We also plot the lines of constant 
ATcou, which is the time left until the universe collapses: H{tQ + ATcou) = 0. We find 
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Figure 3. Confidence levels at 68.3% (light grey inner contour) 95.4% (medium 
grey contour) and 99.73% (dark grey outer contour) are shown in the m- 
00 plane for the potential V{(j)) = Vocos{m(j) / \/\^) . Here m is in units of 
^/SirGVo/S, and 0o is in units of the reduced Planck mass Mp = -^3/871^. 
The best-fit point in each plot is marked as a star. The per degree of 
freedom at the best-fit is 1.050 for r^om = 0.3 and 1.047 for Oom = 0.4. For 
both figures, in the region to the right of the thick solid line, parameter values 
are such that the matter density never reaches the present value, hence this 
region is disallowed by observations. 

Table 2. Time until the end of acceleration, ATcnd- Time until collapse, ATcoii- 







ATend (Gyr) 










(Gyr) 






0.3 


9.8 


34.9 125.7 


698 


.4 


18.2 


55.9 


153.6 


69^ 


iA 


0.4 


9.8 


40.5 111.7 


698 


.4 


18.2 


48.9 


125.7 


69^ 





that the time to collapse is significantly larger than the time to the end of acceleration. 
At the 95.4% confidence level the lower bound to collapse is given by ATcou ~ 18 Gyrs 
for m < 2. 

In Figure El we show the evolution of the deceleration parameter q and the matter 
density Q^a for DDE driven by the cosine potential. We see that models allowed by 
the SNe data show a fairly large spread in the time span during which the universe 
accelerates. Although the time when the universe starts to accelerate is fairly close 
for different parameter values, the future epoch when the universe stops accelerating, 
g(to+ATend) = 0, varies widely between models. The same can be said for the time when 
the universe collapses, which also varies significantly between models. Interestingly, for 




^0 ^0 
Figure 4. A magnified part of the Fig. |31 with (dashed) fines of constant 
^^end (upper panel) and constant ATcou (lower panel) added. In the upper 
panels, the time ATend is measured from the present epoch to when the universe 
stops accelerating: q{to + ATend) = 0. The values of ATond for the dashed 
curves (from top to bottom) are listed in Table I^J (from left to right). For 
both r^om = 0.3 and riom = 0.4, the minimum time taken for the deceleration 
parameter to rise to zero is ATend ^ 0.7H^^ ~ 10 Gyrs (at the 95.4% confidence 
level). For the lower panel, the values of ATcou for the dashed curves (from top 
to bottom) are listed in Table |21 (from left to right). For both r^om = 0.3 and 
^Om = 0.4, the minimum time to collapse is ATcou — 1.3-ff|^^ ~ 18 Gyrs at the 
95.4% confidence level (we assume Hq = 70 km/s/Mpc). In the region to the 
right of the thick solid curve the matter density never reaches its present value 
of Oom = 0.3 (left panel) and fiom = 0.4 (right panel), therefore this region is 
disallowed by observations. 
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Figure 5. The evolution of the deceleration parameter q and the matter density 
ftra is shown for four different DDE models corresponding to different choices of 
m and (pQ in the DDE potential V{(j)) = Vocos(m0/-v/T^) (^^Om = 0.3). Time t is 
in units of y/S/SnGVo. The models have parameter values: m = 1.0, (po = 0.6 
(A), m = 1.0,00 = 0.2 (B), m = 0.74, (/>o = 0.23 (C), m = 1.0, (po = 1.2 (D). 
Models A,B,C are allowed by supernova observations at the 95.4% confidence 
level. The dashed line D in both panels shows the time evolution of q and 
ilin for a DDE model with m = 1.0, (po = 1.2. This model is disallowed by 
observations since the matter density always remains larger than 0.3 (see figure 
The horizontal dotted line in the left panel (g = 0) divides this panel into 
two regions. In the upper region g > and the universe decelerates, whereas 
g < in the lower region in which the universe accelerates. The points of 
intersection of q = with A,B,C show the commencement and end of the 
acceleration epoch in these models. The horizontal dotted line in the right 
panel marks the present epoch when Jlom = 0.3. The solid circles in both left 
and right panels show the epoch when the potential energy of the scalar field 
falls to zero. Note that this occurs after the universe stops accelerating. 

most allowed models, drops to exceedingly small values before the potential becomes 
sufficiently negative to initiate collapse. 

4. Discussion and Conclusions 

In this paper we have examined the possibility that the dark energy responsible for the 
acceleration of the universe decays with time. Models of decaying dark energy have 
recently been proposed in connection with supergravity and M-theory fTli Our 
analysis of decaying dark energy is carried out for two potentials. The first, defined in 
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Eq. (0), is a faithful local representation of a potential in the vicinity of its maximum 
value and therefore describes a very general situation in which the potential responsible 
for dark energy has one (or more) maxima. (In fact Eq. (0) generically describes a dark 
energy potential near its maximum value.) We follow the expansion dynamics of the 
universe as the scalar field rolls down this potential and, for a large region in parameter 
space, compare the resulting value of the luminosity distance with observations of high 
redshift type la supernovae. Our results show that the decaying dark energy model 
is consistent with SNe data at reasonably small values of m and 0o- We see that 
the potential declines to half its maximum value on a time scale not shorter than 
ATi/a ~ O.GHq^ (ATi/2 ^ 8 Gyrs if Hq = 70 km/s/Mpc) at the 95.4% confidence 
level. 

Next, we extend this analysis to the cosine potential (jHl) which holds the possibility 
of becoming negative as the field evolves. The cosine potential has some important new 
features not shared by 'tracker' quintessence potentials. Due to its ability to grow more 
negative with time, the cosine potential allows for the eventual cancellation between 
the positive density of dark matter and the negative density of dark energy, because of 
which the universe collapses at late times. It is important to note that collapse occurs 
generically even in the case of a spatially flat universe, as pointed out in [TT] . Comparing 
our model universe with supernova observations we determine the time, measured from 
the present epoch, to when the universe collapses. At the 95.4% confidence level, the 
time to collapse can never be smaller than ATcou ~ 1.3Hq^ if m < 2 (ATcou ^ 18 Gyrs 
for Ho = 70 km/s/Mpc). 

For both dark energy models (P) & © our current accelerating epoch is a transient 
which ends once the dark energy density falls below a certain critical value. For the 
cosine potential the time until the acceleration of the universe ends {q(to + ATend) = 0) 
is always larger than ATend — 0.7ifJ'^, which is about 10 Gyrs if Hq = 70 km/s/Mpc 
(at the 95.4% CL). Thus a DDE universe will continue accelerating for at least 10 Gyrs 
if i^o = 70 km/s/Mpc and fiom — 0.3. About 8 Gyrs later, this universe will collapse 
and head towards a 'big crunch' singularity. 

Finally, it is worth pointing out that the decaying dark energy models discussed 
in this paper are devoid of horizons, which could be an attractive feature from the 
viewpoint of string/M-theory [Zj. 
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